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We examine a stochastic noise process that has a decohering effect 
on the average evolution of qubits in the quantum register of the solid 
state quantum computer proposed by Kane M. We consider the effects 
of this process on the single qubit operations necessary to perform 
quantum logical gates and derive an expression for the fidelity of these 
gates in this system. We then calculate an upper bound on the level 
of this stochastic noise tolerable in a workable quantum computer. 
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I. INTRODUCTION 



The process of computation by quantum logic, so called quantum computation, 
has recently been shown to be far more powerful for solving certain classes of problem 
than is classical computation |]-[3|]. The superiority stems from the ability of the 
quantum analogue of bits, qubits, to maintain coherence between different classical 
states. This allows the quantum computer to perform computations over many 
classical input states at once, giving a quantum computer an exponential increase 
in speed for solving certain problems. In order to exploit this advantage in quantum 
computation it is vital that the evolution of the qubits be not only coherent, but 
precisely known to the operator. In this paper we examine the evolution of qubits 
under the influence of a stochastic noise process, the effect of which is to make the 
exact evolution of the system uncertain. We thus consider an ensemble of qubits and 
calculate their average evolution. We find that the effect of the noise is to produce 
a decay of the average phase coherence of the qubits in the quantum register and 
to depolarize qubits undergoing single qubit operations. This ensemble decoherence 
manifests itself as a decay of the fidelities of the quantum operations the qubits are 
undergoing. We calculate this fidelity and use it to determine an upper bound on the 
level of stochastic noise that the computer can tolerate yet still operate successfully 
within the limits set by current error correcting codes. 

II. THE KANE SOLID STATE QUANTUM COMPUTER 

Throughout this paper we will be considering a solid state quantum computer 
(QC), as proposed by Kane In this system the qubits are simply spin | 31 P 
nuclei, in a silicon substrate. The system is subject to a background magnetic field 
oriented in the z direction, B z . At low energies the effective Hamiltonian for the 
nucleus-electron system is given by @J^] 
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H B B z a z e - g n /i n B z a* + Aa e .a n . 



(1) 



A = \-KHBgn^n\^{^)\ 



2 and \ip(0) | 2 is the probability density of the electron wave 



function evaluated at the position of the nucleus. We can alter A by applying a 
voltage to a so called "A-gate" situated above the nucleus . This applied volt- 
age shifts the electron wave function away from the nucleus and thus reduces A. 
If A « 1B z \ib the electron spin states are separated in energy by a factor of ap- 
proximately HBl{gn^n) = 1633.8 greater than the nuclear states. Thus the nuclei 
can be manipulated without significantly altering the electron's state. We therefore 
consider as our quantum computing basis a sub-space of the entire Hilbert space 
spanned by the Hamiltonian Eq(|l]), namely the states | { 0), | J, 1). This corre- 
sponds to the basis of nuclear states with the electron in its ground state. The 
effective Hamiltonian in this sub-basis to first order in A/(/j,bB z ), up to a constant, 
is given by 



For convenience we have omitted the n subscript on the Pauli operator as we shall 
do for the rest of the paper. We can tune the Lamour frequency of the nucleus via 
the A-gate bias; 



where V is the applied A-gate voltage, rj = 5n x 10 7 Hz/V and Aq is the value of A 
when V = 0. 

The single qubit operations are implemented by the application of an oscillating 
transverse magnetic field. This field can be brought into resonance by tuning 7 
to satisfy the resonance condition uj = 2B Z ^/Ti. In this way specific qubits can be 
operated on without affecting the rest of the qubits in the register. The Hamiltonian 
for single qubit rotations then becomes H sqr — H + H(t), where 



H = B zl a z . 



(2) 



7 = -9nHB ~ 



A - V V 
B z 



(3) 
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H(t) = -B ac g n fx n (cos(ujt + <f>)a y - sm{u)t + <p)a x ). (4) 

Here the phase factor <fi determines the axis of the rotation. Without loss of gener- 
ality we set = and consider only y rotations. We convert to a frame rotating at 
this resonance frequency by transforming to the interaction picture, and find 

H sqr = -B ac g n fi n a y . (5) 



III. DECOHERENCE OF THE QUANTUM REGISTER 

We now consider the effect of a stochastic white noise in the applied A-gate 
voltage, that is we write the voltage signal 

V(t)=V (l+A(t)), (6) 

where A(t) describes a white noise process |J. We can thus write 

A(t)dt = V\dW(t), (7) 

where dW(t) is the Wiener increment, and \f\ scales the noise. We can calculate 
A by integrating Eq(|) over the duration of the voltage pulse, r, to give the pulse 
area: 

T(r) = V r + AT(r). (8) 

Here AT(r) is a Gaussian random variable with a mean of zero and a variance VqXt. 
The ratio of the rms value of the fluctuations in the pulse area to average pulse area 
is given by: 

We find then, that the Hamiltonian for a qubit in the quantum register, that is a 
qubit not undergoing an operation (B ac = 0), in the presence of this noise is given 
by 



H=B,(<y + £{t))o*. (10) 

Here gives the stochastic fluctuations in the Lamour frequency of the qubit 
caused by the noise, it is related to the noise in the voltage signal by 

?« = fW (id 



We define £(t)dt = ^dW(t), where 

'^ 2 A. (12) 



v B z 

Let us transform into a frame rotating at the Lamour frequency, that is transform 
to the interaction picture as we have already done for the Hamiltonian Eq(|5]). In this 
picture a qubit in a noiseless quantum register does not evolve at all, the Hamiltonian 
is zero. However if we include the white noise, the evolution is generated by the 
Hamiltonian 

H = B z Z(t)° z - (13) 

From this we can form the Ito stochastic differential equation for the density operator 
in the interaction picture 

dffl = ^VtiW(t)[a*,p(t)] 
in 

-^W°*,m]]dt (14) 

Because we do not know the precise history of the evolution, we take an average 
over noise histories and thus calculate the evolution of the average density operator. 
Using (dW{t)) = we find 

d ^ = -^[a\[a\p av {t)\\- (15) 

Master equations with this nested commutator structure have been studied exten- 
sively in the field of quantum optics Let us define polarization vector by 



Pav{t) = |(1 + P(t).a). We then consider the evolution of the average polarization 
vector in the quantum computation basis. In this notation the information about 
coherence is contained in the x and y components of the polarization vector. The 
magnitude of the polarization vector gives the purity of the state; if the state is 
completely pure then \P\ = 1, in a completely mixed state the magnitude is and 
partially pure states have a magnitude between these two extremes. The evolution 
leads to an exponential decay of the x and y components, indicating a loss of phase 
coherence, while the z component remains unchanged: 

dP z (t) 



dt 

dP x , y (t) _ 2eB 2 z 



0, (16) 
P*M (17) 



dt h 

Thus phase coherence is destroyed but population probabilities are conserved. These 
equations have the solution 

P z (t) = P z (0), 

P x Jt) = exph^P^O), (18) 

in which we can explicitly see the exponential decay of phase coherence. We can see 
how this decoherence of the average density operator effects the operation of the QC 
by calculating the average fidelity of the operation. This gives the probability that 
the state evolves as we expect it to if there were no noise, and can be obtained by 
calculating the trace of the product of the evolved density operator and the density 
operator we would expect from a noisless evolution, in this case the zero time density 
operator. 

F = Tr\p av (t)p av (0)}, 

= \(l + PM 2 + (^(0) 2 + P„(0) 2 )exp[^p^]). (19) 

We can see that the fidelity depends on the initial state of the system, specifically 
on how much phase coherence the initial state possessed. If the qubit is initially in a 
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classical state the system cannot decohere and we get a fidelity of 1, a perfect gate. 
Generally however, the initial state will be some kind of superposition of classical 
states and coherence will be destroyed, leading to a loss of fidelity The worst case 
is when the initial state is a maximum superposition, the system decoheres, and the 
fidelity is given by 

F = i(l + expF^]). (20) 

In this case the fidelity eventually decays to a limiting value of |. Because the 
process is only phase destroying, a measurement of a z , that is a projection onto the 
basis of classical states will yield the same result as in the noiseless case. In this 
sense we can say that the classical information has been retained, the qubit has been 
converted to a classical bit. 



IV. FIDELITY OF SINGLE QUBIT OPERATIONS 

To calculate the effect of the stochastic fluctuations in the A-gate voltage bias on 
a qubit undergoing a single particle operation we consider the example of a rotation 
around the y-axis. The inclusion of white noise transforms the Hamiltonian Eq([5|) 
to 

H sqr = i{t)B z <J z - B ac g nf i n a y . (21) 

We can then form the Ito stochastic differential equation 

w,x B ac g n p n -/.ni,, . B z y/edW(t) 
W) = tz — [cr y , p[t)\dt + — [a ,p[t)\ 

-^rl* z ,l* z rp(t)}}- (22) 

Again we average over the noise histories and find 
dp~av{t) B ac g n p n eB z 

dt = ^ W\Pav{t)\ ~ 7^-2 W i W ,Pav[t)\\. (23) 
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This gives a set of coupled differential equations for the polarization vector compo- 
nents of the average density operator 

dP x (t) _ ~2B 2 e 2B ac g n y n 

dt ~ h 2 x[) h nw ' 

dP y {t) _ -2Ble 



dt h 2 yyh 

dP z (t) _ 2B ac g n fi, 

dt h 
that can be solved to give 



Py(t), 

Px{t), (24) 



P x (t) = exp[-g^]{(cosh(g) - ^sinh(|))P,(0) - 2 ^"^\ inh(g)P,(0)}, 
n hah ah 

P y (f)=exph^]P tf (0), 

P,(t) = exp[^]{(cosh(g) + ^sinh(^))P,(0) + ^^sinh(^)P,(0)}. 

(25) 

Here we have defined a = \j B^e 2 — A(B ac g n fi n h) 2 . In this case we see that all the 
components of the polarization vector decay to zero, thus leaving a totally mixed 
state regardless of the initial state. This is known as a depolarizing process, not 
only is the phase coherence lost, but the population probabilities become uniform 
and the qubit is equally likely to be in the classical |0) or |1) state. The action of 
the single qubit gate is to mix the P x and P z components of the polarization vector. 
The dephasing process causes the decay of the P x and P y components and these two 
processes combine to cause the depolarization. The two processes define two time 
scales in the system, the first t ov = p — — is the time it takes to perform a typical 
single qubit logic operation, a Hadamard gate. The time scale of the dephasing 
process is given by r^ ec = ^§p- Obviously a functioning quantum computer requires 
T o P / T dec « 1- To zeroth order in this ratio we find that Eqs(^) give 

P x (t) = exp[^^]{cos( ~ 2 ^^V .(0) + sin(^^^)P 2 (0)}, 
Py(t) = exph^P^O), (26) 



P z (t) = exp[-^]{cos( ^ a ^"> z (0)-sin( 7^ n )P,(0)}. 
We can now calculate the fidelity of these noisy operations 
F = Tr[p av (t)p(t)} 

= \(l + ex P [^^]P y (0r + eM Z §^](PM 2 + PM 2 ), (27) 

where p(t) is the output density operator for a noiseless operation. In this case 
the fidelity decays to | regardless of the input state, however the rate of decay is 
dependent of the initial state. We see that the rate of fidelity loss is not faster 
than is the case for a qubit in the quantum register and in the worst case they 
are equal. The worst case occurs, for y rotations when the qubit is initially in a 
a y eigenstate, this means that the gate cannot rotate the state out of the basis in 
which it decoheres. 



V. TOLERANCE TO NOISE 

To calculate the level of noise in the voltage signal that the quantum computer 
can tolerate in performing a typical single particle operation, the Hadamard gate. 
There is still some debate as to how much error a QC can tolerate and still function 
usefully, even with error correcting codes. The estimates range between an error 
probability of 5 = 1CT 6 — 1CT 4 per qubit per operation 0,0. For the purposes 



of this calculation we will hedge our bets and use the limit 5 = 10~ 5 . We use the 
operating parameters prescribed by Kane H; B z = 2T and B ac = 0.001T and an 
A-gate bias of IV, a value at the top of the bias range. We then find that using the 
worst case fidelity function given by eq(|20|) the error probability is given by 

6 = 1- F 

= ^(l-expP^]), (28) 

" 'dec 



which substituting in our limit for 5 gives 



^<2x 10~ 5 . (29) 

7~dec 



This limit justifies our zeroth order approximation in obtaining eqs (|26| , |2T|) . We now 
require 

e < Bac ^ nh x 4 x 10- 5 . (30) 



Using eq(|l^) we find 



A = — x 4 x 10 > 31 



and so we get a limit on the acceptable noise in the voltage signal. To implement a 
Hadamard gate requires a pulse of duration t = r op , thus we find that the ratio of 
the rms fluctuations in the pulse area to the mean pulse area is restricted by 



F(t H ad) " nr]V h 
which under the operating conditions of the Kane computer gives 

AT(t H ad) 



(32) 



f(t H 



< 1.4 x 10" b . (33) 



ad ) 



VI. CONCLUSIONS 

We have found that the stochastic white noise process which causes dephasing of 
qubits in the quantum register becomes a depolarizing process for qubits undergoing 
rotations in the quantum computer. Under similar operating conditions the fidelity 
of these qubit rotations decays more slowly for certain highly coherent input states, 
than does the fidelity of the register. We find that in order to satisfy currently 
accepted limits of on error probability per qubit operation we must keep the ratio 
of the rms fluctuations in the pulse area to the intended pulse area Ar fa°^"" < 
1.4 x 10~ 6 . 
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